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f^ ' We study in the kr factorisation framework jet production at LHC energies. In par- 

f^ , ticular we are interested in valence quark and gluon initiated jets. The calculation of 

CNJ ' the relevant hard matrix element is presented. A CCFM like evolution equation for 

valence quark distribution is solved and the cross section for valence quark and gluon 
initiated process is calculated using Monte Carlo event generator CASCADE. 

, ■ 1 Introduction 

Collisions at LHC will hopefully answer many so far unanswered questions about Standard 

"^ \ Model and beyond. Here we focus on details of jet production at LHC due to valence quark - 

gluon interaction. This process is of interest since it may shed light on so called phenomenon 
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CIh' of saturation and it may test different factorisation schemes. This process will be measured 
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at central and forward calorimeters (Castor in CMS). 

In first step we calculate hard matrix element for for g*q —^ gq, then we calculate the total 
cross section and show the differential cross sections as functions of Et and ?/ distribution of 
produced jets. To calculate jet production we use CASCADE Monte Carlo event generator 
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■^ ■ 2 Matrix element for g*q — > gq 

CN 

f~^ ■ We begin our study by calculating matrix element for g*q — > gq (this has been calculated in 

^D \ [2] but in a different contest). In our studies we assume Regge kinematics which means that 

we are going to work in regime where the total center of mass energy squared of proton- 
proton collision is much larger than any other scale involved in the collision. This practically 
means: 

s:^ s, \i\, \u\ 



H . where: s = {pA +pb)'^,s ~ {k + qY,u ~ (k — q')'^,i = (fc — fc')^ (see Figure 1). To calculate 

the amplitude we have to sum expressions for diagrams in Figure 1. We proceed by using 
Sudakov decomposition of four momenta: 

k = XgPA + ZgPB + kr 
q^Xg 

k' = Xg:pA + ZgipB + k'j. 

q' = Zq^PA + XqipB + q'r 

Where Xq is longitudinal momentum fraction of incoming quark, Xg is a longitudinal mo- 
mentum fraction of incoming gluon {xg ~ 10~^), fcr is its transverse momentum. Xqi is 
a longitudinal momentum fraction of outgoing quark, q'j, is its transverse momentum. Xgi 
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Figure 1: Diagrams which contribute to QCD Compton scattering in Feynman gauge. In 
calculation we put for proton A an auxiliary quark impact factor and gluon with momentum 
flow k is of shell. Proton B is factored out and the quark with momentum q is considered 
to be on-shell 

is a longitudinal momentum fraction of outgoing gluon, k'j. is its transverse momentum. 
In calculations of matrix element the interesting kinematical regime is where the incoming 
gluon carries low longitudinal momentum fraction {xg ~ 10~^) and we can make the eikonal 
approximation to the current. On the other hand since the incoming quark carries high 
momentum fraction {xq ^ 10""'^) we take the exact expression for the current. The resulting 
formula for the matrix element squared for g*q —^ gq for unpolarized quarks and giuons, 
after summing and averaging over colors of final and initial state particles, reads: 



\Mf = {ATrfa 



X^s'^ix'^ +X'^) f s{8Xq +Xg) -uiSx'g +Xq) 
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From the expression for the matrix element squared we can see that in addition to term 
formally equal to the coUincar result we get a term proportional to k^ . This matrix element 
has interesting structure of singularities. Singularities appear in 4 different points in the 
phase space when one of the final state partons four-momentum is collinear with a four- 
momentum of a quark in the initial state (taking also the quark with the momentum pA 
into account). 

3 Numerical results 

The cross section in kx factorisation approach can be schematically written as: 

da{pp ^gqX)^ Aixg.k'i, ^i^) ® M{g* q^gqX)® Q{xq,ql, ^x^) 

where ® stands for integration over longitudinal, angular and transversal degrees of freedom, 
yui, ^2 are the factorisation scales. In the formula above A(x g , k"^^, ii\) is the unintegrated, 
gluon density obtained by solving the CCFM evolution equation and Q{xq,qf, 112) is the 
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Figure 2: Final quark peudorapidity distribution 
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Figure 3: Final gluon peudorapidity distribution 



unintegrated valence quark density obtained from CCFM-like equation which is needed 
for technical reasons in Monte Carlo simulation. The initial valence quark distribution 
is provided by CTEQ 6.1 set. In our calculations we use running coupling constant and 
we make a lOGeV cut on transverse momenta of out going jets. This cut is motivated 
by the experimental setup where jets at lower momenta are hardly measurable. This cut 
removes the coUincar configuration singularities of the matrix clement. As a first observable 
we calculated total cross section for two jet production in considered kinematical region 
and we obtained 0.l7mb. As other observable wc investigate the rapidity distributions of 
the produced jets. The results of the calculation of the rapidity distributions of outgoing 
partons is presented in Figure 2. As expected, the quark initiated jet is in the forward 
rapidity region while the gluon initiated jet is in the central rapidity region. We make here 
a comparison to PYTHIA Monte Carlo prediction. The difference between distributions 
is due to missing sea quark contribution in CASCADE. The other interesting observable is 
the transverse energy Et (Figure 4) and rj distributions (Figure 5) of produced jets. The 
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Figure 4: Transverse energy distribution of produced jet in Castor calorimeter range. 
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Figure 5: Pseudorapidity distribution of produced jet in Castor calorimeter range. 

first of this distributions is of particular importance because it shows visible difference due 
to different assumptions about underlying physics. For example one sees that CASCADE 
including the off shell matrix element favors a harder spectrum than PYTHIA which is based 
on coUinear factorisation. The 77 distributions generated by CASCADE and PYTHIA differ 
because missing sea quarks in our approach. However, this difference is small in region of 
our interest (large Et)- 
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